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We use the world's data on elastic electron-proton scattering and calculations of two-photon 
exchange effects to extract corrected values of the proton's electric and magnetic form factors over 
the full range of the existing data. Our analysis combines the corrected Rosenbluth cross 
section and polarization transfer data, and is the first extraction of Ge and Gm including explicit 
two-photon exchange corrections and their associated uncertainties. In addition, we examine the 
angular dependence of the corrected cross sections to look for possible nonlinearities of the cross 
section as a function of e. 
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I. INTRODUCTION 

As one of the most fundamental observables which 
characterize the composite nature of the nucleon, elec- 
tromagnetic form factors have over the past few decades 
provided considerable insight into the nucleon's internal 
structure [J 0) S 0] ■ In the commonly used one-photon 
exchange approximation to electron-nucleon scattering, 
the form factors depend only on the four-momentum 
transfer squared, —Q^. At low the Fourier trans- 
forms of the form factors give information on the charge 
and magnetization distributions of the nucleon. 

The expectation from perturbative QCD is that the 
dependence of the Sachs electric, Ge{Q^), and mag- 
netic, Gm{Q^), form factors should be the same at large 
Q^, and early data indeed suggested approximate scal- 
ing of the ratio Ge/Gm with 0]. These data were 
analyzed using the Rosenbluth or longitudinal-transverse 
(LT) separation technique, in which the form factor ra- 
tio is extracted by examining the elastic cross section as 
a function of the scattering angle, 9. The cross section 
at 6* = 180° depends only on the magnetic form factor, 
while the cross section at smaller angles is a combination 
of magnetic and electric contributions. At large Q^, the 
contribution from the electric form factors is small, and 
here the technique has reduced sensitivity to G^;. 

Recent experiments at Jefferson Lab [a, 0, S] , utilizing 
the polarization transfer (PT) technique to measure the 
ratio Ge/Gm, found the surprising result that Ge de- 
creases more rapidly than Gm at large . This indicates 
a substantially different spatial distribution of the charge 
and magnetization of the nucleon |3| , raising questions as 
to the impact of angular momentum and relativistic ef- 
fects on the nucleon structure 10, 11, 12, 13, 14, 15,, 16]. 
The virtue of the PT technique is that it is sensitive only 
to the ratio Ge/Gm, rather than the sum of a large mag- 



netic and small electric contribution, and so does not suf- 
fer from the dramatically reduced sensitivity to of the 
LT separation at large Q^. However, the two techniques 
disagree significantly even in the region where they both 
yield precise results [itJ, and a puzzle has existed about 
the origin of the discrepancy. 

An accurate determination of the proton electric and 
magnetic form factors is also important as these impact 
on our knowledge of other quantities whose extraction is 
sensitive to G^; and Gm- In some cases, uncertainties in 
the proton form factors will be magnified, yielding even 
larger corrections in other observables. 

A number of recent theoretical studies of two-photon 
exchange (TPE) in elastic e» scattering have been per- 
formed [il lii, 0, [IE ii. Ill 111, (see also Ref. HI 
for a recent review). These indicate that TPE effects give 
rise to a strong angular-dependent correction to the elas- 
tic cross section, which can lead to large corrections to the 
LT-extracted Ge/Gm ratio. In fact, the results of quan- 
titative calculations based both on hadronic intermedi- 
ate states I23I l2jl [25l . [27| and on generalized parton 
distributions [20, l21[ , provide strong evidence that TPE 
effects can account for most of the difference between the 
LT and PT data sets. 

While the TPE analyses are very suggestive of the reso- 
lution of the problem, there exists some residual discrep- 
ancy between the TPE-corrected data sets. Moreover, 
the effects of the TPE correction on the extraction of G^; 
and Gm have in practice been estimated using a linear 
approximation to the correction over an arbitrary angular 
range ^19. . ,22i] , which can yield significant uncertainties in 
evaluating the effect of TPE on the form factor extrac- 
tions, especially at high Q^. An accurate determination 
of Ge and Gm requires analyzing the elastic scattering 
data taking into account the TPE effects directly on the 
cross sections. This is the primary aim of the present 
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paper. 

Our approach will be to apply the TPE corrections di- 
rectly to the cross section data, estimate the uncertain- 
ties in the TPE corrections, and compare TPE-corrected 
results to polarization measurements. We then use the 
corrected cross section results, combined with polariza- 
tion transfer measurements, to extract the best possible 
values and uncertainties for Ge and Gm of the proton. 

The outline of this paper is as follows. In Sec.|TT]we re- 
view the formalism of elastic ep scattering in the presence 
of TPE effects and discuss the TPE corrections which 
we apply to the data. We also study the effect of TPE 
on the ratio of electron-proton to positron-proton elas- 
tic cross sections. In Sec. Illli TPE corrections are ap- 
plied to the cross section and polarization transfer mea- 
surements, and a combined analysis of all measurements 
is performed. We provide parameterizations of the cor- 
rected form factors, as well as a direct fit to the elastic 
cross section, without correction for TPE effects. Finally, 
in Sec. IIVI we summarize our results and discuss future 
work. 



II. ELASTIC ep SCATTERING 

A. Born approximation 

In the one-photon exchange or Born approximation the 
differential cross section is given in terms of the Born 
amplitude A^o by: 
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where CMott is the cross section for scattering on a point 
particle, E and E' are initial and final electron ener- 
gies, respectively, r — (5^/4M^, and a — e^/47r is the 
electromagnetic fine structure constant. In performing a 
Rosenbluth separation, it is convenient to work with the 
reduced cross section, an: 
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where e = [l + 2(1 + r) tan^ (^'/2)] is the virtual pho- 
ton polarization parameter. This decomposition allows a 
direct separation of Ge and Gm from the e dependence 
of (Tfl. In the Born approximation the form factors are 
functions of only a single variable, Q^. Effects beyond 
one-photon exchange introduce additional e dependence. 

The polarization transfer method involves the scatter- 
ing of a polarized electron beam from an unpolarized tar- 
get, with measurement of the polarization of the recoiling 
proton. In the Born approximation the ratio of the trans- 
verse to longitudinal recoil polarizations yields [28j : 
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where is the proton magnetic moment, and Pt {Pl) is 
the polarization of the recoil proton transverse (longitu- 
dinal) to the proton momentum in the scattering plane. 
It should be noted that both Pt and Pl individually 
depend only on the ratio Ge/Gm, with Pt (Pl) be- 
ing strongly (weakly) dependent on Ge/Gm- In most 
experiments, however, the ratio of the two polarization 
components is used, with Pt providing sensitivity to the 
form factor ratio, and the ratio to Pl yielding a result 
that is independent of the overall beam and target po- 
larization and/or analyzing power of the detector. One 
has similar sensitivity to Ge/Gm in measurements of the 
beam-target asymmetry [2^ [s^, HH . 



B. TPE: Unpolarized ep scattering 

With the inclusion of radiative corrections to order a, 
the elastic scattering cross section in Eq. ([1]) is modified 
according to: 
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where 5^q represents one-loop corrections, including vac- 
uum polarization, electron and proton vertex, and two- 
photon exchange corrections, in addition to inelastic 
bremsstrahlung for real photon emission ^3^. 

As discussed in Refs. 19, 22], the amplitude for the 
one-loop virtual corrections, Aii, can be written as the 
sum of a factorizable term, proportional to the Born am- 
plitude A^Oi and a non- factorizable component, Aii- 
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The ratio of the full, 0{a) cross section to the Born cross 
section can then be written as: 
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with the correction given by: 
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The contributions to the functions f{Q'^,s) from the 
electron vertex, vacuum polarization, and proton vertex 
terms depend only on Q^, and hence do not affect the 
LT separation, aside from an overall normalization fac- 
tor. Of the factorizable terms, only the IR-divergent TPE 
correction contributes to the e dependence of the virtual 
photon corrections p^ . 

The non- factorizable terms contained in A^i arise from 
the finite nucleon vertex and TPE corrections, which de- 
pend explicitly on nucleon structure. For the proton ver- 
tex correction, the e dependence is weak, and will not 
significantly affect the LT analysis For the inelastic 
bremsstrahlung cross section, the amplitude for real pho- 
ton emission can also be written in the form of Eq. ^ . 
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In the soft photon approximation, in which one of the 
two exchanged photons is taken to be on-shell, the full 
amplitude is completely factorizable. A significant e de- 
pendence arises due to the frame dependence of the an- 
gular distribution of the emitted photon. These correc- 
tions, together with external bremsstrahlung, contain the 
main e dependence of the radiative corrections, and are 
accounted for in standard experimental analyses. They 
are generally well understood, and in fact enter differ- 
ently depending on whether the electron or proton are 
detected in the final state [13, [11] . 

The only remaining term at 0{a) that can introduce a 
non-negligible e dependence, namely the non-factorizable 
part of the TPE contribution, is typically not accounted 
for in cross section extractions. Its effects can be included 
by considering the interference of the total TPE and Born 
amplitudes, according to: 

2^{mI M2-,] 
— L L tQ\ 



where the TPE amplitude M.2-y in principle includes all 
possible (off-shell) hadronic intermediate states. 

In typical experimental analyses of electromagnetic 
form factor data, radiative corrections are implemented 
using a formalism based on the Mo & Tsai (MT) prescrip- 
tion |^[3^. While the detailed implementation has been 
improved over the years d, [13], the treatment of TPE 
has remained unchanged from the original work. The 
TPE effects are partially included in the MT approach 
by approximating the TPE amplitude by its infrared (IR) 
divergent part. 

A more detailed examination of the loop integrals [ssj 
yields the same IR-divergent contribution as Mo & 
Tsai [111, 5ir(MT), but also yields an IR-finite contri- 
bution. The IR-finite part, which is usually neglected in 
the standard data analyses, was found in Refs. |19l.[2^ to 
have a significant e dependence, and is explicitly included 
here. To isolate the effect of the additional TPE contri- 
butions on the data, Blunden et al. [l^ HI] consider the 
difference: 

A = - SmiMT) , (9) 

in which the IR divergences cancel. We will follow this 
approach and therefore reference the e dependence of the 
full calculation of 62-y with that of 5ijt(MT). 

The results for the difference, A(e,(5^), between the 
full calculation and the MT approximation are most sig- 
nificant at low £, and essentially vanish at large e. At 
the lower values, A is approximately linear in e, but 
significant deviations from linearity are observed with in- 
creasing Q^, especially at smaller e (2^ . 

The effect of TPE on the LT extractions of G^; is great- 
est at large since the cross section here is not very sen- 
sitive to Ge- However, the TPE corrections to the cross 
section have a relatively weak dependence, and are 
therefore still significant at low values, where many 



high-precision measurements of nucleon or nuclear struc- 
ture are performed that require accurate knowledge of 
the proton form factors. 



C. TPE: Theoretical uncertainty 

For the dominant nucleon elastic TPE contribution, 
the main uncertainty arises from the input form factors 
at the internal 7*7ViV vertices in the loop diagrams. Be- 
cause the intermediate state nucleon is off-shell, in prin- 
ciple the "half off-shell" form factors here need not be 
the same as the free nucleon form factors, and can explic- 
itly depend on the intermediate nucleon four-momentum. 
The off-shell dependence is of course unknown, but in 
practice it is sufficient, at least at low Q^, to approxi- 
mate these by the free form factors. 

For the calculation in Ref. [22| it was necessary for 
technical reasons to parameterize the internal form fac- 
tors by sums of monopolc functions, which were fitted 
to realistic parameterizations of form factor data. In the 
actual analysis [22l ] the form factor parameterization was 
taken from the global fit in Ref. [3g]. For comparison, 
more recent fits to data [1^, [13] were also used, and the 
sensitivity of the TPE correction A(e, Q^) to the particu- 
lar form was found to be negligible up to ^ 10 GeV'^. 

In the present analysis we use the more recent global 
form factor parameterization from Ref. [39] with Ge 
constrained by the PT data. A three-monopole fit to 
this parameterization provides a very good fit for Ge 
and Gm over the range where both the form factors are 
constrained by data. Above = 6 GeV^, the three- 
monopole fit to Ge becomes unconstrained, leading to 
unstable results for the TPE corrections. Here we in- 
stead use a dipole approximation for the form factors, 
which guarantees a smooth extrapolation to high Q^, 
and matches onto the asymptotic high-Q^ behavior ex- 
pected from perturbative QCD. The difference between 
the TPE results at = 6 GeV^ with the "realistic" and 
dipole form factors is not significant, which is consistent 
with the findings in Ref. |22l ]. Similar calculations from 
Ref. [1^ find results that are in very good agreement with 
those of Ref. 

Estimates of the contributions to A(e, Q^) from higher- 
mass intermediate states have been made by a number 
of authors [H, [13, [iH [i^. [H, [13]. In a recent analy- 
sis, Kondratyuk et al. [23| evaluated the contribution of 
the A(1232) resonance, which is known to play an im- 
portant role in hadronic structure, between = 1 and 
6 GeV'^. The A(1232) intermediate state contribution 
was found to be smaller in magnitude than the nucleon 
contribution, with an opposite sign at backward scatter- 
ing angles where the TPE effects are largest. For realistic 
choices of the 7iVA vertex, and couplings obtained from 
a hadronic model analysis of Compton scattering off the 
nucleon [i^, the magnitude of the A(1232) contribution 
was found to be around 1/4 of that of the nucleon at 
low and intermediate . The A(1232) therefore tends 
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to cancel some of the TPE effect from the nucleon elastic 
intermediate state, modifying the cross section by be- 
tween — 1% and 2%. At larger the magnitude of the 
A(1232) contribution increases, along with that of the 
nucleon elastic, although above ~ 10 GeV^ the relia- 
bility of the calculation is more questionable. 

Effects of higher mass, spin- and isospin-1/2 and 3/2 
resonances were also estimated recently by Kondratyuk & 
Blunden [l^ . The excited states included in that analysis 
were the Pn iV(1440) Roper resonance, the D13 iV(1520), 
and the odd-parity S\\ iV(1535) for isospin-1/2, and the 
S'31 A(1620) and D33 A(1700) isospin-3/2 resonances. 
The photocouplings for these states were taken from the 
model [4^ of nucleon Compton scattering at energies up 
to the first and second resonance regions, and a universal 
dipole form factor was applied for each of the 7-/V 
resonance transitions. In addition, the assumption was 
made that the transitions are mostly magnetic, with the 
electric being much smaller, and the Coulomb couplings 
taken to be zero. 

The overall scale of the higher mass contributions is 
approximately an order of magnitude smaller than the 
nucleon and A(1232) for ;^ 6 GeV^. Moreover, these 
alternate in sign, with the Pn, and D33 contributions 
being mostly negative (as the nucleon), while the and 
S'31 mostly positive. For larger this model, as with all 
models based on hadronic degrees of freedom, is unlikely 
to be reliable. 

Clearly the inclusion of the higher mass hadronic in- 
termediate states is considerably more model dependent 
than the nucleon elastic or even the A(1232) contribu- 
tion, in view of uncertainties in the various photocou- 
plings, the transition form factors, and the efhcacy of 
representing high energy off-shell nucleon excitations by 
states of zero width. Effects of overlapping resonances, 
as well as the multihadron continuum, or non-resonant 
background, are expected to be increasingly important 
at larger . One could in principle constrain the high 
mass spectrum phenomenologically by Compton scatter- 
ing data at high . Unfortunately these are not yet 
available, and one must rely on theoretical guidance. 

From the point of view of quark-hadron duality, one 
expects that at large the intermediate state spec- 
trum can be saturated either by including a large set of 
hadronic resonances, or by summing over quark degrees 
of freedom (46j . Although the approaches are in principle 
complementary, if is large enough a quark level cal- 
culation may provide a more efficient description. Such 
calculations have been performed in the context of gen- 
eralized parton distributions (GPDs) 20, .21] within the 
"handbag" approximation, in which the scattering of the 
virtual photon takes place off a single quark in the proton. 
This should be the dominant (leading twist) contribution 
at high , although at intermediate Q^, processes where 
the virtual photon interacts with different quarks in the 
proton (higher twist) are still important. 

One of the uncertainties in the calculation is the choice 
of GPD, and in Ref. [2ll a Gaussian valence ansatz is 



used. The calculation has a small but noticeable differ- 
ence when using their modified Regge GPD, and has a 
similar dependence on the assumed quark mass. In ad- 
dition, for the leading twist approximation to be appli- 
cable, both and the invariant mass squared s of the 
intermediate state are required to be large. In practice 
the constraint Q^,s > is used, which restricts the 
kinematic reach of the calculation to moderate and large 
e. For Q2 _ 2-6 GeV^ one is restricted to e > 0.2-0.4. 
The results show significant e dependence with clear de- 
viations from linearity at intermediate and large e, and 
a weak dependence. The overall trend of the e de- 
pendence is similar to that in Ref. [2^ . but smaller in 
magnitude. At high there appears little e dependence 
for e > 0.5, which yields little modification to the LT ex- 
traction of Ge/Gm dominated by data at large e. At 
the highest points the partonic contribution explains 
about 1/2 of the discrepancy between the LT and PT 
results. 

The indications from both the hadronic and partonic 
higher mass calculations are therefore that there is an 
additional effect at the level of a few percent on top of 
the nucleon elastic contribution, albeit with a sizable the- 
oretical uncertainty. Given the difhculty in obtaining a 
reliable quantitative estimate of the higher mass contri- 
butions to TPE at high , rather than rely on a specific 
model, we take a more phenomenological approach. 

At low Q^, the correction is approximately linear with 
e. While the calculations suggest increasing nonlinear- 
ities, they differ in the detailed e dependence, and the 
data are consistent with very small nonlinearities p^ . 
We assume therefore a linear e dependence for the form 
of the extra TPE correction. We make an estimate of the 
dependence of the additional terms from the calcula- 
tions of the higher resonance states and the GPD- 
based model [2l|, and apply this as the additional TPE 
correction. For < 1 GeV^, we do not apply any ad- 
ditional TPE effect, as all indications are that it is very 
small. For > 1 GeV^, we take a dependence that 
grows slowly with Q^: 
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{Q^ > 1 GeV^), (10) 



with in GeV^, and apply a 100% uncertainty. This 
yields a correction linear in s that decreases the e = 
cross section by 1% at = 2.2 GeV^ by 2% at 
4.8 GeV^, and by 3% at 10.6 GeV^. This correction is 
then added to the calculated 52-y 

As a further check, we compare the extraction of 
Ge/Gm from the TPE-corrected cross section data and 
the polarization transfer measurements, which is dis- 
cussed in Sec. IIII Bl and shown in the bottom panel of 
Fig. [21 below. The additional correction brings the high- 
LT Ge/Gm points into good agreement with the PT 
data, although the error on the corrected LT data is con- 
siderably larger at high Q^. 

While the effect of the TPE corrections on Gm is 
smaller, it can potentially have a large impact on global 
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fitting. For instance, the TPE corrections for Ge shift 
around 10 - 15 data points by 1 — 2cr, but for Gm there 
are of order 30 - 40 points affected by a 2 — 3cr shift. 
In addition, even though the TPE corrections to Gm are 
less dramatic than those for Gs, the uncertainties in the 
TPE effects when extracting Gm at large will still be 
important. 



D. TPE: e~p vs. e+p 

While TPE corrections modify the unpolarized cross 
section, it is difhcult to isolate TPE corrections from the 
Born cross section experimentally without independent 
knowledge of the proton form factors. Only deviations 
from the linear e dependence predicted by Eq. ([2]) can 
be observed, and in practice these are found to be very 
small [475. However, direct experimental evidence for the 
contribution of TPE can be obtained by examining the 
ratio of e^p and e~p cross sections. The Born amplitude 
changes sign under the interchange e~ e"*", while the 
TPE amplitude does not. The interference therefore has 
the opposite sign for electron and positron scattering, 
which will show up in the experimental ratio = 
<^e+p/^e-p ~ (1 — 2A), where A is defined in Eq. 
Interference between electron and proton bremsstrahlung 
also yields a small difference, but this is corrected for in 
the extraction of R'^^'^ . 

Although the current data on elastic e~p and e^p scat- 
tering are sparse, there are some experimental constraints 
from early data taken at SLAG |48l, ;49|] jCorneU ^50], 
DESY [53 and Orsay [H (see also Ref. ^ and refer- 
ences therein). The data are predominantly at low 
and at forward scattering angles, corresponding to large 
e [e ^ 0.7), where the calculated 27 exchange contribu- 
tion is small 1%). Nevertheless, the overall trend in 
the data reveals a small enhancement in R"^^*^ at the 
lower e values [sl]. 

Figure [1] shows the extracted ratio, i?^^^ , compared 
with the calculation of Ref. [l^l (filled diamonds), and 
with the additional high mass contribution in Eq. (jlOp 
(open diamonds). The calculation is in good agree- 
ment with the data, although the errors on the data 
points are quite large. Clearly better quality data at 
backward angles and moderate Q^, where an enhance- 
ment of up to ~10% is predicted, would be needed for 
a more definitive test of the TPE mechanism. There is 
a planned experiment to perform a precise (~1%) com- 
parison of e~p and e^p scattering at = 1.6 GeV^ 
and £ « 0.4 at the VEPP-3 storage ring [11] ■ An ex- 
periment to provide broader e and coverage [1^ at 
Jefferson Lab is approved to make such measurements 
up to = 2 — 3 GeV^ using a beam of e+e" pairs 
produced from a secondary photon beam. 
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FIG. 1: (Color Online) Ratio of elastic e^p to e~p cross sec- 
tions. The data (crosses) are taken from Ref. [53] and refer- 
ences therein. The filled diamonds show the ratio using the 
TPE calculation of Ref. 22], while the open diamonds in- 
clude the additional contribution in Eq. (|10p . which modifies 
the high-Q^ correction at low e. 



E. TPE: Polarization observables 



Naively, one expects the corrections to the spin- 
dependent cross sections to be of the same order of mag- 
nitude (<; 5%) as the corrections to the unpolarized cross 
sections. For unpolarized scattering, this yields a large 
correction to the extracted value of G^; at large Q^, where 
the total contribution from Gb is also small. The polar- 
ization measurements are directly sensitive to the ratio 
IJ-pG E I G M ^ and so the TPE corrections to this ratio will 
be of the same magnitude as the correction to the cross 
section; the large magnification of the effect observed in 
the Rosenbluth separations does not occur for polariza- 
tion measurements. 

Both the hadronic [13] and partonic [2l[ approaches 
can be used to estimate the TPE corrections to the re- 
coil polarization measurements. While the calculations 
are not strictly valid in the same kinematical regimes, 
both predict the corrections to be extremely small, ^ 1%, 
except for large and low e values. The corrections to 
the existing data are less than 2% for all settings ex- 
cept for the largest point, where the corrections are 
of order 5%, which is still much smaller than the sta- 
tistical uncertainty of the measurement. However, the 
sign of the correction is different in the two approaches. 
As the sign of this correction is not known empirically, 
and all indications are that it is very small compared to 
the experimental uncertainties, we do not apply a TPE 
correction to the polarization measurements. An experi- 
ment 'sg'I at Jefferson Lab will map out the e dependence 
of the polarization transfer at fixed Q^, which will pro- 
vide the first direct information on the TPE correction 
to polarization transfer. 
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III. GLOBAL ANALYSIS 

In this section we perform a combined analysis of cross 
section and polarization data, including both the stan- 
dard radiative corrections and the TPE contributions de- 
scribed above. In Ref. [13, the form factors were ex- 
tracted after applying Coulomb distortion corrections, 
but only cross section data were included, and the anal- 
ysis was limited to very low values. This is the first 
global analysis of form factors in which TPE effects are 
included in a consistent way from the outset. We apply 
the TPE corrections to the data, and repeat the global 
analysis of Ref. [39|], with a few important modifications: 

• We include additional low cross section data ISSl . 
[5i,[60|, more recent cross section results [s^ l6ll |. 
and the full set of polarization transfer and beam- 
target asymmetry measurements 

• We use TPE corrections following the formalism of 
Blunden et al. [l^l, rather than the phenomeno- 
logical correction used in Ref. [s^. At high 
(above 1 GeV^) we supplement this with an ad- 
ditional TPE correction (Eq. with a 100% 
uncertainty, to estimate the contributions from the 
high mass intermediate states and related (conser- 
vative) uncertainties. 

• We use a different fitting function, based on Ref. , 
and include high-Q^ constraint points in the fit. 
The goal is to have a fit valid at both very low 
and the highest values of the existing data, with 
a sensible extrapolation to higher values. 

In addition to extracting parameterizations of and 
Gm from the global fit, we also take the corrected cross 
section data in small bins, and perform direct LT sep- 
arations. This allows a test of the consistency between 
the TPE-corrected cross section results and the polar- 
ization measurements. Using the same bins we also 
perform a combined fit to the e dependence of the cross 
sections and the polarization data, to extract both values 
and uncertainties for Ge and Gm- 

Finally, we provide a parameterization of the TPE- 
uncorrected cross section measurements. This parame- 
terization is nearly identical to taking the TPE-corrected 
form factors and applying the TPE calculation used here, 
but does not require an explicit calculation of the TPE 
effects. 



A. Global fitting procedure and results 

The combined analysis uses the same general approach 
outlined in Refs. [13, HI]. As in the previous analyses, 
we remove the small angle data from Ref. [^, and ap- 
ply updated radiative correction factors for lepton loop 
diagrams to some of the older experiments which did 



not include these corrections. In Ref. [35| a high preci- 
sion measurement of the elastic cross section was aimed 
specifically at the ratio (IpGe/Gm- Their quoted sys- 
tematic uncertainty is separated into three components: 
normalization, point-to-point, and a linear e dependent 
correction which impacts ^IpGe/Gm, but does not yield 
random fluctuations or nonlinearity relative to the lin- 
ear fit. We use only normalization and point-to-point 
uncertainties, and to include these data in the fit we in- 
creased the point-to-point systematic uncertainty from 
0.45% to 0.5-0.55% (depending on Q^) such that the un- 
certainty on ij-pGe/Gm using the inflated point-to-point 
uncertainty matched the uncertainty as extracted in the 
more detailed analysis. 

After applying the TPE correction to the raw cross 
sections, we perform a global fit to all of the cross sec- 
tion, polarization transfer, and beam-target asymmetry 
measurements using the form: 

GE,GM/fip = ' ■ ■ (11) 

We take n = 3, which yields 8 fit parameters for each 
form factor, along with a normalization factor for each 
independent set of cross section measurements. The fit 
includes 569 cross section and 54 polarization transfer 
data points, and yields a reduced of 0.770. Taking 
n = 2 yields similar fits, with a reduced of 0.802, 
and changes that are below the uncertainties in the form 
factors; Gm changes by ^1% up to = 10 GeV^, Ge 
by less than 3% up to = 5 GeV^. The normalization 
factors are consistent between the n = 2 and n = 3 fits, 
with an RMS variation below 0.2%. The polarization 
transfer data yield a contribution to x^ of 52 for the 54 
data points, indicating good consistency with the cross 
section data. The parameters from the combined fit are 
given in Table HI and shown as solid lines in Fig. [31 The 
parameterization is valid up to at least = 6 GeV^ for 
Gb, and Q2 ^ 30 GeV^ for Gm- 



TABLE L Fit parameters for the extracted proton electric and 
magnetic form factors, using the parametrization of Eq. (1111) . 



Parameter 


Gm/ fip 


Ge 


ai 


-1.465 


3.439 


a2 


1.260 


-1.602 


as 


0.262 


0.068 


bi 


9.627 


15.055 


62 


0.000 


48.061 


63 


0.000 


99.304 


64 


11.179 


0.012 


65 


13.245 


8.650 



The fit function in Eq. ITTI) is chosen as it gives rea- 
sonable behavior in both the limits of low and high 
Q^. We constrained the parameters in the denominator 
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{bi in Eq. (fTTj) ) to be positive, to avoid fits where both 
the numerator and denominator pass through zero at the 
same Q^, yielding narrow divergences in the fit. Finally, 
extra high "data" points were added to the fit to pre- 
vent uncontrolled behavior in the large limit. The 
high constraint points were Gm = 0.7 Go, where 
Gd = [1 + QV(0.71 GeV2)]-2 is the dipole form factor, 
at Q2 = 50,100,200, and 400 GeV^, and = at 
= 10, 15, 20, and 25 GeV^. In both cases, the uncer- 
tainties were taken to be 50%, 100%, 150%, and 200% 
of Gd, respectively. This has a negligible effect in the 
range of the data, but prevents "extreme" behavior in 
the region where the data do not provide a meaningful 
separation between Ge and Gm- 



B. Rosenbluth analysis 

The combined analysis yields a fit for Ge and Gm, as 
well as normalization factors for the various cross section 
measurements. However, it does not give a clear indi- 
cation of the consistency between the cross section and 
polarization data. In this section, we use the normaliza- 
tion factors from the combined global fit to normalize the 
individual cross section data sets, and perform a global 
analysis using only the cross section data binned in Q^. 
The results from the TPE-corrected cross section analysis 
can then be compared directly to the PT measurements, 
to see to what extent the TPE corrections resolve the 
discrepancy. 

We begin by taking cross section data in small bins 
and interpolating each cross section to the mean value 
using the fit from Tab.|Tl The cross section is known well 
enough that the uncertainty associated with the inter- 
polation is small. The bin size is chosen to keep the 
interpolation correction small, typically at the 1% level, 
and an additional uncertainty equal to 3% of the correc- 
tion is applied to the corrected point. We then perform 
an LT separation to extract Ge/Gm- We determine a 
systematic uncertainty associated with the scale factors 
by taking each data set and varying its normalization by 
the estimated scale uncertainty. Note that because we 
have used the global fit to estimate the normalization 
factors, we have a smaller scale uncertainty than quoted 
for the individual measurements. We take the scale un- 
certainty after the global fit to be the smaller of 1% or 
half of the originally quoted normalization uncertainty. 
For the analysis including TPE corrections, we also take 
the additional TPE correction arising from the high mass 
intermediate states and vary it by ±100% to estimate the 
additional TPE uncertainty for large . 

In Fig. [2] we show a comparison of the PT data (filled 
diamonds) and the LT separation data (open circles). 
The data sets are completely independent, except for the 
use of the normalization factors from the global fit of 
Sec. nil Al The top panel shows the global LT analy- 
sis without applying TPE corrections, while the bottom 
is with TPE corrections. The uncertainties in the LT 




[GeV^] 



FIG. 2; (Color Online) Ratio ^pGe/Gm extracted from po- 
larization transfer (filled diamonds) and LT measurements 
(open circles). The top (bottom) figure shows LT separations 
without (with) TPE corrections applied to the cross sections. 



separation increase significantly at large values when 
the TPE corrections are applied. This is in part due to 
the uncertainty assigned to the high TPE corrections, 
but mainly due to the fact that the LT separation is a 
measure of [Ge/GmY rather than Ge/Gm, so a reduc- 
tion in the slope yields a constant absolute uncertainty 
in [Ge/GmY, which maps into a larger absolute uncer- 
tainty in Ge/Gm- 

Note that for very low Q^, the TPE correction does not 
go to zero, as was also pointed out in Ref. [6§|. In fact, the 
correction to the unpolarized cross section becomes zero 
for « 0.3 GeV^, and changes sign for lower values. 
This appears to be a largely model-independent result, 
which persists even for the case of point-like nucleons. 
Thus, even the extremely low extractions are modified 
by TPE. With the TPE correction applied, the average 
value of ^pGb/Ga/ for < 0.2 GeV^ goes from 0.988(4) 
to 0.997(4), showing that Ge and Gj\/ have consistent 
low-Q^ behavior. 

One of the unique features of the TPE correction is 
that it introduces nonlinearity in the e dependence of the 
cross section. Although part of the TPE contribution is 
linear in s, observation of nonlinearity in the data would 
provide direct evidence of TPE effects in elastic scatter- 
ing. To quantify the amount of nonlinearity in the data. 
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we fit the reduced cross section to a quadratic in 
Ref. j47l: 

an = Po[l+ Pi{e - 0.5) + P2{e - Q.hf] , 



where P2 represents the fractional e-curvature parame- 
ter, relative to the average {e = 0.5) reduced cross sec- 
tion. With the inclusion of TPE corrections, the average 
nonlinearity parameter, (^2), is found to increase from 
1.9 ± 2.7% to 4.3 ± 2.8%. While the extracted nonhn- 
earity increases with the TPE corrections, it is not large 
enough to be considered inconsistent with P2 = 0. In 
addition, the results from Ref. [47j are dominated by 
higher points, where we do not include nonlinearities 
in the TPE contributions from higher mass intermediate 
states. Including the single-experiment LT separations 
from the new low data sets used in this analysis, we 
find (P2) = 2.8 ± 2.4% (after TPE), still generally con- 
sistent with no nonlinearities. 



C. Extraction of Gb and Gm from global analysis 

In this section we extract individual Ge and Gm points 
and uncertainties over the full range where the form 
factors can be separated. The analysis follows that of 
the corrected cross section data in the previous section, 
but now we include the PT measurements in each bin 
as part of the fit. The results for < 6 GeV^, where 
Ge and Gm can be separated, are given in Table HIl and 
shown in Fig. [3l 

For the combined fit, we fit Ge and Gm to the com- 
bination of cross section and polarization transfer data 
in each small bin. For the cross sections, we use the 
TPE-corrected cross section measurements, and normal- 
ize each data set using the scale factors found in the 
global fit (Sec. IIII After making the initial fit for Ge 
and Ga/, we scale each data set by the estimated nor- 
malization uncertainty (as in Sec. IIIIBp to find its con- 
tribution to the systematic uncertainty, and add these 
in quadrature to determine the total uncertainty in Ge, 
Gm and the ratio due to the normalization uncertainties. 

In addition to improving the overall precision, combin- 
ing the cross section and PT results has the added benefit 
of decreasing the correlation in the uncertainties in G^; 
and Gm- The Rosenbluth separation tends to yield a 
large anti-correlation between the uncertainties for Ge 
and Gm, and thus an enhanced uncertainty on the ratio. 
The PT data measure the ratio directly, thus dramati- 
cally reducing this correlation. Therefore, in Tab. |TT1 we 
provide values and uncertainties for both the individual 
form factors and the form factor ratio. 



D. Extraction of Gm at high 

In the extraction of Gm for > 6 GeV^, the value 
of Ge is not known, so that an additional assumption 
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FIG. 3: (Color Online) Extracted values of Ge and Gm from 
the global analyses. The open circles are the results of the 
combined analysis of the cross section data and polarization 
measurements fSec. IIII"C] Tab. The magenta crosses are 
the extracted values of Gm (Tab. IIIIll for the high re- 
gion, where Ge cannot be extracted. The solid lines are 
the fits to TPE-corrected cross section and polarization data 
(Sec. IIII A|) . The dotted curves show the results of taking 
Ge and Gm to be ^/ol and y/or, respectively, from a fit to 
the TPE-uncorrected reduced cross section ( Appendix K)) . i.e. 
the value one would obtain using only cross section data and 
ignoring TPE. 



is required in order to extract Gm- We extract Gm 
under two different assumptions for the ratio Ge/Gm- 
First, we assume that the Ge term is negligible above 
6 GeV^, which would be the case if G^ approached zero 
and then stayed small. Second, we assume a linear fall- 
off, UpGe/Gm = 1-0.135 (g2-0.24), from Ref. 0. Up 
to « 14 GeV^ this yields a smaller contribution from 
Ge than in previous analyses, where it was assumed that 
IJ-pGe/Gm — 1. 

At higher Q^, the linear fit yields \ijLpGe/Gm\ > 1, and 
thus a larger Ge contribution, almost 10 times what was 
assumed in the inital analysis of the Sill, et al. data [T^I 
at = 30 GeV^ . This yields a significant change in the 
dependence of the reduced cross section. Instead of 
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the electric contribution decreasing from 6% to almost 
zero as one went from 5 to 30 GeV^, the linear decrease 
in implies almost no contribution at 5-9 GeV^, and 
a maximum contribution of « 6% at the highest Q^." 

The final quoted values for Gm are taken as the av- 
erage of the values obtained assuming \IpGeIGm = 0, 
and iMpGe/Gm = 1 - 0.135(Q2 _ 0.24), with half of the 
difference taken as the associated uncertainty. Note that 
the additional TPE component meant to estimate the 
high mass intermediate state contributions, as discussed 
in Sec. Ill C[ yields a significant correction (and uncer- 
tainty) of 3-4.5% to the cross section over the range of 
the high data. Table IIIII shows the extracted val- 
ues of Gm based on the above assumptions, along with 
the values of Gb needed to reproduce the TPE-corrected 
cross sections when combined with the quoted values of 
Gj\/. Because we average the extracted Gm values, and 
the cross section correction due to Gb depends on G\, 
this corresponds to taking the average value of G\ rather 
than the average Ge- Above (5^=10 GeV^, the only data 
included is from Sill, et al. [701. There are also data from 
an earlier SLAG experiment [71[ , but both the scale and 
statistical uncertainties are much larger than for the later 
measurements, and so these data are not included. 



IV. CONCLUSION 

We have performed the first global analysis of elastic 
electron-proton scattering data taking into account two- 
photon exchange contributions and their associated un- 
certainties. The analysis combines the corrected Rosen- 
bluth cross section and polarization transfer data, and 
the corrected form factors Ge and Gm have been ex- 
tracted over the full range of of existing data, up to 
g2 « 6 GeV2 for Ge and 39 GeV^ for Gm- 

The TPE corrections applied are based on the hadronic 
model of Ref. [23| for the nucleon elastic intermediate 
state, but with improved input form factors at the in- 
ternal vertices. These have been supplemented by in- 
elastic contributions estimated from recent calculations 
with explicit excited N* intermediate states |23, 27, 41], 
and from GPD-based partonic calculations [20] . The un- 
certainty in describing the (off-shell) intermediate states 
is the main source of model dependence in the analysis. 
This uncertainty is more tractable, and relatively mild, 
at low , but becomes larger at high values. For the 
higher-mass intermediate state contributions, which are 
more important for the high data, we assign a 100% 
uncertainty to the estimated inelastic TPE correction. 

The resulting TPE corrections to Gb are significant in 
the region 2 ^ ^ 6 GeV^, where the LT and PT data 
are in most striking disagreement, and here they bring 
the LT data into good agreement with the PT results 
for Ge/Gm- The corrections to Gm are smaller, but 
can be a few percent at moderately large values. We 
provide a convenient parameterization of the corrected 
Ge and Gm form factors, as well as of the reduced elas- 



tic cross section, parameterized by effective form factors 
which contain TPE effects. 

While several recent quantitative studies have demon- 
strated that much of the disagreement between the LT 
and PT results could be accounted for by TPE, none 
of these calculations can be said to be complete, and 
none fully resolve the discrepancy at all values. The 
hadronic calculations used here have the least model de- 
pendence at lower Q^, but require assumptions about 
the additional high corrections which have large un- 
certainties. There could in principle also be larger than 
expected TPE corrections from high mass states to the 
PT experiments, although the indications from all exist- 
ing studies suggest that these are likely to have a small 
effect on the form factor ratio. Additional measurements 
of TPE corrections will help to refine these calculations 
by directly determining the TPE effects to both the cross 
section and polarization measurements. 

One such measurement is experiment E04-116 (ssj . us- 
ing a beam of e^e~ pairs produced from a secondary pho- 
ton beam at Jefferson Lab, which will make simultaneous 
measurements of e^p and e~p elastic cross sections up to 
^ 2 GcV^. The e^p/e^p ratio is directly sensitive 
to TPE effects. A proposal to perform a precise (~ 1%) 
comparison of e~p and e+p scattering at = 1.6 GeV^ 
and e « 0.4 has also been made at the VEPP-3 storage 
ring [5I]. 

With the inclusion of two-photon exchange, we also 
expect nonlinearities in the e dependence of the cross 
section. With the existing data, the results do not show 
a significant nonlinearity before or after applying TPE 
corrections, and so do not yet set significant limits. A 
recently completed Jefferson Lab experiment [72i] will 
provide an accurate measurement of the e dependence 
of the elastic ep cross section, with sufficient sensitivity 
to test the calculated nonlinearities. Additional upcom- 
ing experiments (sgI . pTSj will extend polarization transfer 
measurements to higher values, and will examine in 
detail their e dependence. The increase in will allow 
for a clean extraction of Ge and improved uncertainties 
on Gm up to ~9 GeV^, while the e dependence will be 
the first measurement sensitive to TPE in polarization 
transfer. 

On the theoretical front, the largest uncertainty arises 
from poor knowledge of the contributions to the TPE 
amplitude from high mass intermediate states, beyond 
the nucleon elastic contribution. This introduces signif- 
icant model dependence in the correction, especially at 
large Q^. It can be reduced by using phenomenological 
input for the virtual Gompton scattering amplitude at 
intermediate and high values. In future it would be 
desirable to merge the low hadronic calculations with 
the high partonic approach, and develop a frame- 
work which enables TPE corrections to be consistently 
described over the entire range of data. The present 
analysis should at the very least serve to highlight this 
need. Finally, we may expect in future more informa- 
tion on virtual Gompton scattering to come from lattice 
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QCD, although calculation of four-point functions of this 
type is currently still in its infancy. 



APPENDIX A: PARAMETERIZATION OF 
ELASTIC CROSS SECTION 

In this appendix we provide a convenient parameter- 
ization of elastic cross section. One can use the fits to 
Ge and Gm in Sec. IIII CI to calculate the Born cross 
section, and then apply radiative corrections, including 
TPE, to predict the measured cross section. However, 
this requires an explicit calculation of the TPE correc- 
tion to the cross section, consistent with the calculation 
used to correct the cross section measurements. It is use- 
ful for some applications to work directly with the elas- 
tic cross section, without necessarily requiring knowledge 
of the fundamental form factors. We therefore provide 
a separate parameterization of the elastic cross section 
and uncertainties without TPE corrections applied to the 
data. This will provide a simple and reliable, model- 
independent parameterization for the cross section with- 
out direct reference to TPE effects. 

We parameterize the full elastic reduced cross section 
(Born + TPE corrections) by: 

^Born+TPE ^ ^^2 (q2^ ^) ^ sF^iQ^s) , (Al) 

where Fm and Fe are effective magnetic and electric form 
factors, respectively, which absorb the effects of multiple 
photon exchange. Note that this is a purely phenomeno- 



logical fit form, rather than a true representation of the 
effect of two-photon exchange on the form factors, which 
introduces a new term p^ . In the Born approximation 
these obviously approach the usual Sachs form factors: 

FUQ\e) Gm(Q') , F,{Q^,e) -> G£(Q2) . (A2) 

In principle, the form factors i^e and F,„ depend on 
both and £, but as calculations and existing data [i^l 
indicate extremely small deviations from linearity, we can 
at present safely neglect this additional e dependence, 
and use the form from Eq. pT|) when fitting Fm and Ff.- 
We repeat the global fit described in Sec. IIII A( using 
only uncorrected cross section data. We keep the nor- 
malization constants from the different data sets fixed to 
the values obtained from the global fit, as these represent 
out best estimation of the true normalization factors of 
the cross section measurements. The results of the fit 
for the cross section in Eq. (jAip in terms of the effective 
form factors and Fm are given in Table HVl 
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TABLE II: Extracted values for the corrected form factors 
relative to the dipole form, and the ratio [x^GeIGm- 



'^[GeV^] 






D 


G 


M /iHpGo) 


(j-pGe/Gm 





.007 


1, 


.000±0, 


.006 


0, 


.715±0, 


.393 


1, 


.400±0, 


.712 





.012 


0, 


.996±0, 


.006 


0, 


.966±0, 


.152 


1, 


.032±0 


.154 


0, 


.017 


0, 


.995±0 


.003 


0, 


.984±0 


.031 


1, 


.011±0 


.033 


0, 


.022 


0, 


.993±0, 


.003 


0, 


.953±0, 


.019 


1, 


.042±0, 


.022 


0, 


.030 


0. 


.988±0, 


.007 


0. 


.983±0, 


.052 


1. 


.005±0, 


.059 


0, 


.038 


0, 


,987±0 


.004 


0, 


,985±0 


.013 


1, 


,002±0 


.016 


0, 


.041 


0, 


,974±0 


.007 


1, 


,016±0 


.027 


0, 


.959±0 


.031 





.048 


0, 


.981±0 


.005 


0, 


.982±0 


.009 


0, 


.999±0 


.012 





.057 


0. 


.980±0, 


.005 


0, 


.992±0, 


.010 


0, 


.988±0, 


.012 





.061 


0, 


.990±0 


.008 


0, 


.978±0, 


.015 


1, 


.013±0 


.021 


0, 


.069 


0, 


.977±0, 


.007 


0, 


.991±0 


.016 


0, 


.986±0, 


.021 





.081 


0, 


.977±0, 


.005 


0, 


.980±0, 


.008 


0, 


.997±0, 


.011 





.098 


0, 


.973±0 


.007 


0, 


.977±0, 


.007 


0, 


.996±0 


.010 


0, 


.115 


0, 


.974±0, 


.011 


0, 


.978±0 


.007 


0, 


.996±0 


.016 


0, 


.138 


0, 


.964±0, 


.012 


0, 


.981±0, 


.008 


0, 


.982±0, 


.018 


0, 


.171 


0, 


.986±0, 


.012 


0, 


.977±0, 


.007 


1, 


.009±0, 


.015 


0, 


.199 


0, 


,979±0 


.011 


0, 


.981±0 


.006 


0, 


.998±0 


.013 


0, 


.234 


0, 


.978±0, 


.012 


0, 


.979±0, 


.012 


0, 


.999±0, 


.015 





.273 


0, 


.955±0, 


.010 


0, 


.972±0, 


.008 


0, 


.983±0, 


.014 


0, 


.304 


0. 


.960±0, 


.010 


0. 


.981±0, 


.007 


0. 


.978±0, 


.013 





.350 


0, 


.939±0 


.014 


0, 


.985±0, 


.013 


0, 


.953±0 


.009 


0, 


.390 


0, 


.965±0, 


.011 


0, 


.993±0, 


.008 


0, 


.972±0, 


.017 





.428 


0, 


.961±0 


.015 


1, 


.003±0, 


.014 


0, 


.958±0, 


.010 





.473 


0, 


.970±0 


.011 


0, 


.995±0 


.008 


0, 


.976±0 


.011 


0, 


.528 


0, 


.984±0, 


.013 


0, 


.996±0, 


.009 


0, 


,988±0 


.016 





.584 


0, 


.967±0, 


.013 


1, 


.007±0, 


.007 


0, 


.960±0, 


.016 





.622 


0, 


.969±0 


.014 


1, 


.007±0, 


.007 


0, 


.962±0, 


.020 


0, 


.689 


0, 


.981±0, 


.023 


1, 


.017±0 


.010 


0, 


.965±0, 


.031 


0, 


.779 


0, 


.965±0, 


.011 


1, 


.021±0, 


.005 


0, 


.945±0, 


.013 


0, 


.853 


0. 


.916±0, 


.022 


1. 


.045±0, 


.007 


0, 


.876±0, 


.025 


0, 


.979 


0, 


.933±0, 


.049 


1, 


.045±0 


.017 


0, 


,893±0 


.060 


1, 


.020 


0, 


.920±0 


.017 


1, 


.054±0 


.006 


0, 


,873±0 


.018 


1, 


.170 


0, 


.922±0 


.014 


1, 


.059±0, 


.005 


0, 


.871±0 


.016 


1, 


.330 


0. 


.936±0, 


.047 


1, 


.057±0, 


.010 


0, 


.885±0, 


.051 


1, 


.520 


0, 


.889±0, 


.022 


1, 


.074±0, 


.006 


0, 


.828±0, 


.023 


1 


.740 


0, 


.844±0 


.020 


1, 


.077±0, 


.004 


0, 


.784±0, 


.020 


1, 


.830 


0, 


.856±0, 


.038 


1, 


.076±0, 


.008 


0, 


.796±0, 


.038 


2, 


.070 


0, 


.837±0, 


.038 


1, 


.075±0 


.006 


0, 


.778±0, 


.039 


2, 


.500 


0, 


.744±0, 


.038 


1, 


.077±0 


.005 


0, 


.691±0, 


.037 


2, 


.690 


0, 


.723±0, 


.075 


1, 


.076±0, 


.009 


0, 


.672±0, 


.075 


2, 


.900 


0, 


.676±0, 


.039 


1, 


.079±0, 


.007 


0, 


.626±0, 


.037 


3 


.090 


0, 


.673±0, 


.118 


1, 


.075±0, 


.010 


0, 


.626±0, 


.116 


3 


.240 


0, 


.673±0, 


.112 


1, 


.074±0, 


.010 


0, 


.627±0, 


.110 


3 


.710 


0, 


.652±0, 


.040 


1, 


.068±0, 


.005 


0, 


.610±0, 


.039 


4, 


.110 


0. 


.546±0, 


.059 


1. 


.063±0, 


.007 


0. 


.514±0, 


.057 


5, 


.010 


0. 


.409±0, 


.057 


1. 


.046±0, 


.005 


0. 


.391±0, 


.055 


5, 


.850 


0. 


.287±0, 


.095 


1. 


.027±0, 


.007 


0. 


.280±0, 


.093 



TABLE III: Extracted values for the high Gm results, 
along with the Ge values assumed in the extraction. 



Q [GcV J 


f-rM / [flpUD ) 


Ge/Gd 


7.081 


0.996±0.014 


+0.054 


8.294 


Q.940±0.Q18 


-0.058 


9.848 


Q.908±0.Q13 


-0.191 


11.990 


0.877±0.Q21 


-0.364 


15.720 


0.816±0.028 


-0.629 


19.470 


0.721±0.Q33 


-0.814 


23.240 


0.713±0.Q39 


-1.062 


26.990 


0.689±0.Q47 


-1.273 


31.200 


0.696±0.067 


-1.564 



TABLE IV: Parameters for the fit to the TPE- uncorrected 
cross section of Eq. (|A2|I . using the parameterization of 
Eq. jn]) for Frn and Fe. 



Parameter 


Fm 1 \-^'p 




ai 


-2.151 


-1.651 


£12 


4.261 


1.287 


0-3 


0.159 


-0.185 


bi 


8.647 


9.531 


b2 


0.001 


0.591 


bs 


5.245 


0.000 


bi 


82.817 


0.000 


bs 


14.191 


4.994 



